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The excitation of the axial quasi-normal modes of a relativistic star by scattered particles is 
studied by evolving the time dependent perturbation equations. This work is the first step towards 
the understanding of more complicated perturbative processes, like the capture or the scattering of 
particles by rotating stars. In addition, it may serve as a test for the results of the full nonlinear 
evolution of binary systems. 
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I. INTRODUCTION 



The gravitational interaction of massive bodies in general relativity is associated to the emission of gravitational 
signals, the features of which depend on the nature of the interacting bodies and on the characteristics of their orbits. 
To be studied in detail, very energetic processes like the coalescence of binary systems need the complex machinery of 
the fully non linear theory of gravity coupled to hydrodynamics; this challenging task has been undertaken by several 
scientific groups in the world ||], Q, Waiting for the outcome of these simulations, approximation methods 

which allow to get preliminary information on these processes are extremely useful. An example of the way in which 
approximate methods shed light on much more complex phenomena, is the capture of test particles by a black hole 
(see ref. § for a review). On the assumption that the infalling mass is much smaller than the black hole mass, so 
that its effect can be treated as a perturbation of the black hole geometry, it was shown that the gravitational signal 
emitted in this process exhibits a ringing tail due to the excitation of the quasi normal modes of the perturbed black 
hole. Years later, a fully relativistic numerical simulation of the gravitational collapse of a massive star to a black 
hole showed that the signal emitted in this process also exhibits an exponentially damped sinusoidal tail. The 
experience matured in black hole perturbations suggested that this tail appears because the newly born black hole 
radiates its residual mechanical energy oscillating in its quasi-normal modes. 

A further element of interest in the use of perturbation theory is that the results will serve as a test for fully 
nonlinear calculations 0. For instance, the nonlinear evolution of the head on collision of two black holes |^ has 
proved to be in excellent agreement with the results obtained via a perturbative approach. 

The gravitational signals emitted by massive particles orbiting around rotating and non-rotating black holes have 
been extensively investigated Q, and more recently this study has been extended to the signals emitted when small 
masses in open or circular orbit, interact with a compact non rotating star [^- jT^ . By a direct integration of the 
equations describing the perturbation of the star, it has been shown that both its fluid modes and the spacetime 
modes can in principle be excited, if the scattered mass is allowed to get sufficiently close to the compact star. These 
results have been obtained by Fourier transforming the perturbation equations, and solving the associated boundary 
value problem, with the condition that the solution is regular at the center of the star, and reduces to a pure outgoing 
wave at radial infinity. 

This approach is very powerful when masses are in open or closed orbits around a non rotating star; however, it 
cannot be applied if the mass falls onto the star, due to the lack of information on the interaction between the particle 
and the stellar fluid after the impact. To investigate this kind of problems, a time evolution approach seems more 
appropriate than the analysis in the frequency domain, since it allows to compute the emitted signal at least up to 
when the particle touches the surface of the star. It appears promising also to describe the perturbations of rotating 
stars excited by small masses. Indeed, in that case the coupling between the axial and polar perturbations which 
arises as a consequence of the dragging of inertial frames ( , ) , leads to a set of equations that are not simply 
reducible to wave equations as in the case of a non rotating star. Moreover, the time evolution approach has recently 
been used in the literature to study the scattering of wave packets by non rotating relativistic stars, ( 0, 0, [Q, 
Ipof), where it has been shown to correctly reproduce the spectra of axial and polar perturbations, in a perturbative 
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study of colliding neutron stars |20| (close limit), and to provide proper initial data for stellar perturbations pl| ]. 

These considerations motivate the work done in this paper: we develop a scheme for the numerical evolution of stellar 
perturbations excited by a scattered mass. In order to allow a direct comparison with previously obtained results, 
we integrate the equations describing the axial perturbations of a homogeneous star with increasing compactness, for 
different values of the orbital parameters of the scattered mass, and we choose the same stellar and orbital parameters 
as in ref. where the axial equations were integrated in the frequency domain. 

It should be mentioned that the time-evolution equations do not suffer of the annoying problem of the divergence of 
the source term when the particle transits through the turning point, which causes some difficulties in the integration 
of the Regge- Wheeler and of the Zerilli equations in the frequency domain. 

The plan of the paper is the following. In section II we shall briefly review the time-dependent version of the 
equations of the axial perturbations. In section III we shall present the results of the numerical integration. Our 
energy spectra will be compared with those found in and and the differences will be discussed. In the Appendix 
we shall show that the low-frequency part of the spectrum can be traced back to the bremsstrahlung radiation emitted 
by the scattered mass. 



II. THE EQUATIONS DESCRIBING THE PERTURBED SPACETIME 

In order to describe the non axisymmetric perturbations of a star induced by a scattered particle, we expand the 
perturbed metric in tensor spherical harmonics, and choose the Regge- Wheeler gauge. The perturbed line element 
has the form: 

ds^ = e^^dt^ - e^^^dr^ - r^d-d^ - sin^ Mip^ 

.dYirn r,o 



2sin. 
sin dip 



[hlJt,r)dtdip + h]^{t,r)drdip] 
[h°,^it,r)dtdd + hl^{t,r)drdd] 



(2.1) 



We consider the perturbations of a non rotating star with uniform energy density e =const [g^ . It should be reminded 
that homogeneous stars can exist only if their radius R exceeds 9/8 times the Schwarzschild radius, or R/M > 2.25 
(Buchdahl limit). The source of the perturbed Einstein equations is given by the stress-energy tensor of the particle 
moving along a geodesic of the Schwarzschild spacetime {T{t),R{t), 0(i), $(t)) : 



dTjt) dz^^ dz- 5{r - R{t)) 



(2.2) 



With this source, and assuming that the particle moves on the equatorial plane Q = ^, the axial perturbations can 
be shown to be described by the following equation ( |2^ , ) 



d^Zim{t,r) d^Zir,r{t,r) 



dt^ 



drl 



+ Veir)Zi.rait,r) 



im 



(2.3) 



where 



Zern{t,r) 



-hi{t,r) 



(2.4) 



is a gauge invariant quantity, is the tortoise coordinate and 

Vi = (e^^/r^) [£{£ + l)r + r^{e^p)- 6m{r)] 



(2.5) 



is the potential, which reduces to the Regge- Wheeler potential outside the star, where pressure and density vanish 
(p = e = 0). The two functions Dim and Qem, related to the particle's motion, have the following form: 



Qlrn{r,t) 



i^i.e-'^'' dR 



£2 e"^" 

^—5{r-R(t))-YL^ 



e=f ,<f(t) 



£{£ + 1) r3 dt 



Sir - Rit)) ■ YL^ 



0=f ,*(*) 



(2.6) 
(2.7) 
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where fi is the mass of the particle. 

The one-sided energy spectrum of gravitational waves at infinity will be calculated from the following relation 

dj^i^E E {i-m^+m+mucor, (2.8) 

{=2 rn=~l 

where Zim.{co) is the Fourier transform of the wavefunction Z£m{t,r), evaluated at radial infinity. For the axial 
perturbation, we have to deal only with multipoles with m ~ I — 1,1 — 3, . . . since the multipoles with m — 1,1 — 2, . . . 
correspond to polar perturbations. For the purpose of this paper we will present results only for I — 2 and m — 1. 



III. NUMERICAL METHOD AND RESULTS. 



The numerical procedure we follow to calculate the energy spectrum can be described by the following steps 

• For a specific stellar model and for a given set of values of the energy and of the angular momentum of the 
particle, (E, L^), we calculate the geodesies i.e. the position of the particle, for a discrete number of time steps 
from r oo wp to r — Tt {rt is the turning point of the given open trajectory). 



For the specific geodesic, we estimate the source term of (2.3) for each time and spatial grid point. Since there 



is no back reaction this calculation is done only once at the beginning. 



According to equations (2.6)-(2.7) the position of the particle at each time step is represented by a delta function. 



that we approximate by a narrow Gaussian function 

5(r - « ^e-^'^'-^W)' (3.1) 



/TT 



(this approximation has been also used in |1^ ). The code has been tested for various values of a in order to 
check the way in which the width of the Gaussian affects the results. 



Finally, eq. (2.3) is evolved using a second order finite-difference scheme. The results have been tested for 



various grid sizes in order to achieve the required accuracy. 

In Figure |l| we plot two waveforms corresponding to a particle with orbital parameters (E^L^) = (2.38/i, 12/iM) 
scattered by two homogenous stars, one with compactness R/M = 2.26 and the other with R/M = 2.30. The initial 
part of the signal is similar in both cases, since it is due mainly to the accelerated particle, whereas the subsequent 
oscillations have different frequencies and damping times since they are due to the excitation of the axial modes of the 
two stars. The energy spectrum of the above signals, calculated via eq. (2.8), is shown in Figure ^ On the same figure, 



we mark the location of the real part of the frequency of the axial quasi-normal modes of the two stars, as computed 
in by solving the homogeneous equation as a boundary value problem. It comes out that the lowest trapped 
modes, that have much longer damping times, do not seem to be excited, whereas the higest w-modes are clearly 
excited in both cases. In Figure || we show the energy spectrum for a less energetic particle (E, L^) — (1.01^, i.5fiM); 
here one can easily identify the peaks corresponding to the excitation of the w-modes, but the bulk of the emitted 
energy is due to the quadrupole emission of the accelerated particle. This continues to be true if the central star is less 
relativistic. These results suggest that in order to have a significant excitation of the stellar axial modes one needs 
both the periastron of the scattering orbit to be as close as possible to the stellar surface, and the star to be of high 
compactness. The above features agree with the results of ref. where the axial and polar emission for particles 
scattered by a compact polytropic star was studied in the frequency domain, and with those presented in |13| , where 
the polar mode excitation by scattered particles has been studied with a time-evolution approach. 

Our results also qualitatively agree with those of ref. , as a comparison of our figure 2 and 3 with their figures 
5, 6 or 8 clearly show. However, there is a difference of a factor ~ 27r in the energy spectra, which is probably due to 
a different normalization. More relevant is the difference we find in the energy spectra at low frequency: we observe 
some peaks that do not appear either in the spectra of ref. , or in those presented in ref. . Since this same 
feature also appeared in previous works ( |l0|] and |13|), we have investigated the low frequency behaviour in some 



more detail. We have computed the energy spectrum emitted by the same scattered particles in a semirelativistic 
approximation, by assuming that the particles move along a geodesic of the background geometry, but radiate as if 
they were in fiat spacetime, according to the quadrupole formula. The energy spectra estimated by this procedure 
are shown in Figure | for the star with R/M = 2.26 and for {E, L^) = (2.38/Li, 12^M) and [E, L,) = (1.01^, 4,.b^iM) 
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(details of this calculation are given in the Appendix). We find low frequency peaks which correspond to those in 
figures ^ and ^; this indicates that the low frequency part of the energy spectrum is "shaped" by the bremsstrahlung 
radiation emitted by the accelerated particle. It should be mentioned that one should not expect an exact coincidence 
between the results of the semirelativistic approximation and of the relativistic approach, especially when the orbit of 
the scattered particle is very relativistic. However, the two approaches must tend to the same limit when the particle 
moves on less relativistic orbit, and we have checked that this is indeed the case. 

As a concluding remak, we would like to point out that the calculations presented in this paper prove that the 
numerical evolution of the time dependent perturbation equations is a fast and accurate way to study scattering 
processes occurring in the vicinity of a neutron stars; they allow to obtain detailed information either on the amount 
of energy emitted at low frequency, which can be traced back to the quadrupole emission, and on the excitation of 
the neutron star quasi- normal modes. In a future work we plan to extend the present investigation to the study of 
the excitations of the stellar modes of rotating stars. There, the coupling of the axial and polar perturbations makes 
the whole process more complicated and the numerical evolution of the perturbation equations appears a promising 
approach. 



Helpful discussions with E.Berti and J.Ruoff have been greatly appreciated. We are also grateful for a grant from 
the INFN which has enabled our collaboration. 



In order to understand what is the contribution of the gravitational radiation emitted by the accelerated particle 
to the total energy spectrum, we have used a semi-relativistic approximation introduced in 1971 [ p4| , which assumes 
that the particle moves along a geodesic of the curved spacetime, but radiates as if it were in flat spacetime. On this 
assumption, using the quadrupole formula, we have computed the TT-components of the gravitational wave emitted 
by the scattered mass. As usual, the reduced quadrupole moment is 



where the two-dimensional vector X is the position of the particle along its trajectory in the equatorial plane 
X = cos sin and R{t) and ^{t) are given by the geodesic equations. The expressions of the 

second time derivative of the components of Qki in terms of R{t) and <I>(i) are 
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APPENDIX: THE QUADRUPOLE EMISSION OF THE SCATTERED MASS 




(3.2) 




2R{t)^ + 2R{t)R{t) - 4i?(t)2$(t)2 cos ^{t)^ 



&R{t)^{t)R{t) + 2R{tf^t) sin $(t) cos $(t) 




(3.3) 




vv 



2R{tf + 2R{t)R{t) - AR{tf^{t) cos<^(tf 



+ 



m{t)^{t)R(t) + 2R{tf^{t) sin $(t) cos $(t) 




(3.4) 




xy — 



8R{t)^{t)R{t) + 2R{tf^{t) cos$(t)2 



+ 



2R{tfxa + 2R{t)R{t) - m{tf^tf sin $(t) cos $(t) 



R{tf<l>{t) - AR{t)^t)R{t) 



(3.5) 




(3.6) 



From these expressions the non vanishing TT-components of the emitted wave can be computed as follows 
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hf;^it, x) 



(3.7) 



where Pn 



UiUk is the projector onto the 2-sphere r ~ const, and n is the radial unit vector. The exphcit 



expression of the two independent components is 



rh. 



TT 



Qxx ~ Qyy ) COS^ If + 2Qxy sin (f COS ip (l + cos^ z9) 



Qyy Q z 



cos^ ■& + Q^^ - 



-2cos?9 



[Qxx ~ Qyy^ sin ipCOSip + Q^y (1-2 COS^ ifi) 



(3.8) 
(3.9) 



where i),ip are the polar angles. Using the geodesic equations the Q ajt ) give n by eqs. (3.3)-(3.6) can be computed. 
Then, by Fourier-transforming the metric components given in eqs. ( |3.8| )-( |3^ ) h'^J{uj,r,'d,(p), h'^J{uj,r,d,ip) are 
easily evaluated. Since the energy per unit frequency and unit solid angle is 



dE 



dQdaj IGtt^ 



(3.10) 



where the frequency co is restricted to be positive, the energy spectrum can be computed by integrating over the 
solid angle. It should be pointed out that our convention on the Fourier transform is /(w) = f{t)e'^'^*'dt. 

We have numerically evaluated the energy spectrum for the two sets of values of the orbital parameters considered 
in this paper, i.e. {L^ = 12M, E — 2.38) and (L^ — 4.5Af, E = 1.01) , and the results are plotted in figure |4[ 
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FIG. 1. The waveforms of the signal emitted by a particle with {E, Lz) = (2.38/i, 12fi,M) scattered by a star with R/M = 2.26 
(upper panel) and R/M = 2.30 (lower panel). The excitation of the w-modes is clearly manifested in the signal. 
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FIG. 2. The energy spectrum of the signal emitted by a particle with (E,Lz) ~ (2.38/i, 12/iM) scattered by a star with 
R/M = 2.26 (upper panel) and R/M — 2.30 (lower panel). The sharp peaks correspond to the excitation of the w-modes, while 
the broader peaks at low frequency are mainly contributed by the quadrupole radiation of the particle due to its orbital motion. 
The diamonds correspond to the real part of the known quasi-normal mode frequencies. It seems that the lowest frequency 
trapped modes of the first model of star are difficult to be excited. 



6 




O.O O.I 0.2 0.3 0.4 0.5 O.e 

IS/1 CO 



FIG. 3. The energy spectrum of the signal emitted by a particle with {E,L^) = (1.01/x, 4.5/iM), scattered by a star with 
R/M — 2.26 and with R/Al — 2.30. The sharp peaks correspond to the excitation of the various w-modes of the star, whereas 
the dominant contribution to the emitted energy is that of the quadrupole emission of the scattered particle, which produces 
the broader peaks at lower frequencies. 
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FIG. 4. The energy spectrum computed by the semirelativistic approximation described in the Appendix, is plotted for two 
particles with {E, L^) — (2.38p, 12^M), and {E, Lz) = (1.01/x, 4.5/iM) scattered by a star with R/M — 2.26. A comparison with 
the energy spectra given in figures |and|, indicates that the low frequency emission is contributed by the the bremsstrahlung 
radiation emitted by the accelerated particle. 
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